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A magnetic torque method is proposed that probes the warping and mass gap of Dirac cone surface 
states in topological insulators like 612X3 (X=Se,Te). A rotating field parallel to the surface induces 
a paramagnetic moment in the helical surface states for nonzero warping. It is non-coUinear with 
the applied field and therefore produces torque oscillations as function of the field angle which are 
a direct signature of the surface states. The torque dependence on field strength and angle, the 
chemical potential and the Dirac cone parameters like warping strength and mass gap is calculated. 
It is shown that the latter leads to a symmetry reduction in the fourfold torque oscillations. 
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I. INTRODUCTION 

The strong topological insulators of the 
612X3 (X=Se,Te) family, including the various doped 
systems are characterized by surface states comprising 
a single Dirac cone at the T point of the Brillouin zone 
(BZ)2. These states are protected against the effect of 
disorder, e.g., impurity doping as long as time reversal 
invariance is preservecP"Sl. In the simplest effective 
models the surface states have a helical spin polarization 
orthogonal to the wave vector on the circular Fermi 
surface (FS) which is cut out from the cone at the 
position of the chemical potential. However, it has 
become clear that this picture is oversimplified. When 
the chemical potential is moved away from the Dirac 
point the cone becomes warped and the FS changes from 
circular to hexagonal snowflake shape. This effect has 
been directly observed in photoemissiorf^ and explained 
as the result of higher order terms in the effective surface 
HamiltoniarPIS]. Furthermore ferromagnetic order of Fe 
dopants may lead to the opening of a mass gap in the 
Dirac cone due to breaking of time reversal invariancd^. 

In experiments it is not always easy to discrimina te the 
contribution of surface states from those of the buHp^HH. 
In this work we propose a straightforward and sensitive 
magnetic torque method to observe the physical effect of 
surface states and signatures of their cone structure. The 
schematic setup is shown in the inset of Fig. [2] A mag- 
netic field B is applied and rotated parallel to the surface 
of a disk-like sample. The c-axis is oriented perpendicular 
to the plane. When the chemical potential is away from 
the bottom and top of bulk bands a paramagnetic mo- 
ment due to polarization of helical spins of surface states 
may be induced. Note that in this geometry there is no 
diamagnetic effect, i.e., no Landau quantization of Dirac 
cones. However the spin polarization leads to a moment 
M which is non-collinear to B, giving rise to a magnetic 
torque — V^slM x B = rgz where Vsl is the total 
volume of the surface sheets on both sides of the sample. 
The torque may be measured via piezoelectric elements 
fixed to the disk (inset in Fig. [2]) . In fact recently this 



method has been successfully applied to measure the tiny 
torque oscillations from a collective hidden order param- 
eter in a low carrier heavy fermion systerrfSl. A theoret- 
ical analysis of this effect was given in Ref. [TS] In the 
present case the torque oscillation is entirely due to the 
virtual particle-hole continuum excitations of the warped 
Dirac cone. We propose this method as a promising way 
to investigate the surface states. Furthermore the torque 
amplitude is sensitive to the chemical potential and the 
opening of a mass gap at the Dirac point caused by mag- 
netic surface doping (e.g. by ¥^^. When the moments 
are perpendicular to the surface there will also be an in- 
plane torque (_L z) originating directly from the magnetic 
surface layer. However, it is independent of the direction 
of B and also independent of the chemical potential. 
Because it is perpendicular to the torque from Dirac sur- 
face states it should be possible to discriminate it from 
the latter by an experimental setup as shown in Fig. [2] 

In Sec. |lT] we derive the expression for the surface 
torque based on an effective low energy model for surface 
states. The discussion of numerical results is presented 
in Sec. |III| and the conclusion is given in Sec. |IV[ 



II. MAGNETIC TORQUE AND 
PARAMAGNETIC SUSCEPTIBILITY IN THE 
WARPED DIRAC CONE MODEL 

There are two contributions to the total out-of plane 
torque originating from the bulk {tb) and topological 
surface states {ts) respectively. The latter is obtained 
from the paramagnetic susceptibility tensor x'^ of surface 
states by M = X'^B. Then we can express the surface 
torque contribution aJ^ 

^s(0) = ^l/sLS'[Ax''(0)sin20-2xf^(0)cos20] (1) 

with Ax"^ = Xxx ~ Xyy and Vsl = "^.ASs {A— disk area, 
(5s=surface state localisation length which is of the or- 
der of a quintuple layer thickness in 612X3). Here (f> 
is the azimuthal field angle , i.e., B = B {cos (p, sine/)). 
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FIG. 1: Contour plots of warped and gapped Dirac cones in an 
in-plane magnetic field (0 = 37r/4). To be compatible with 
Eqs. (3]4l we use the convention of kx,ky corresponding to 
TK, TM directions in the hexagonal surface BZ respectively. 
Here intrinsic units {kx,ky) = {bkx,bky) (Table [l| are used 
with m = m/E"; B = 'ysB/E* where the energy scale is E* = 
vp/b — {v%/\)^. In clockwise direction (a-d): {m,B) = 
(0,0); (0.3,0); (0,0.1); (0.3,0.1). The mass term (chosen large 
for clarity) reduces the rotational symmetry while the field 
shifts the Dirac or extremal point of the dispersion. 



The question arises whether can be separated from 
the bulk part. The latter is due to the occupied valence 
band states. If the chemical potential varies within the 
bulk band gap, then tb will be constant whereas ts will 
strongly depend on the chemical potential as discussed 
in Sec. |III| Therefore if the difference oi t = tb + ts 
at chemical potentials fi and —fj,c (inset of Fig. |3| is 
considered, then the surface part may be extracted by 
TsilJ') = t{ii) — t{—^c) since Ts{~fJ.c)—0- The relative 
magnitude of ts to tb cannot be calculated in the present 
effective low energy surface state model used below. A 
rough estimate may be obtained from the following con- 
sideration: If we take a thin-film like sample composed 
of n quintuple layers of 612X3 and assume that the sur- 
face states are localised within one unit one may estimate 
ts/tb ^ 2/(n — 2) for ~ 0. The surface state Dirac 
cones exist for n > 5 (Ref. I16p , then in the optimum case 
one has ts/tb — 2/3. Even for n ~ 10^ ts would still 
be 2% of Tb- Since torques can be measured with rela- 
tive accuracy of 10"'* (Ref. Ill]) there seems to be no real 
obstacle to extract the surface torque contribution. Nat- 
urally one should try to observe this effect in thin films 
as function of the layer number similar to the photoemis- 
sion experiments for samples with n < 50 (Ref. il6|) and 
also to thin film transport measurement^^. 

To identify the effect of surface states on the torque 
we now have to calculate the field angular dependence 
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FIG. 2: Top: Variation of torque amplitudes A\,i for /i = 
0.2, B = 0.1 at (j) — n/8 with (dimensionless) Dirac cone 
warping parameter A = Xq^/fj,c, here m = 0. Bottom: Total 
torque ts for rh — (full line) and m = 0.1 (dashed line). 
For finite m the symmetry under sign change of A is lost. 
Inset shows schematic geometry with disk-like sample fixed 
by piezoelectric elements (grey circles) to measure the torque 
Ts{(j>) under rotating field B. 



of the susceptibility tensor caused by the surface states. 
An effective low energy surface model may be obtained by 
starting from four atomic spin-orbit and crystalline elec- 
tric field split orbitals of quintuple layers in 612X3 with 
appropriate surface boundary conditions^. The form of 
the model HamiltoniarP is dictated by inversion symme- 
try I , time reversal T and the threefold rotation C3. 
Explicitly one obtains: 

+ma^ - js{BxO-x + ByCTy) (2) 

where k± — k^ ^ iky. The first parabolic term i?ok = 
k^/2m* does not contribute to Ax'^ or Xxy ^'^'^ there- 
fore irrelevant for the torque as is shown below. The sec- 
ond term describes the perfect Dirac cone and the third 
term is the third order warping part whose strength is 
given by the constant A. Similar fifth order terms will 
be neglected. The fourth term is due to T reversal sym- 
metry breaking caused e.g. by ferromagnetically ordered 
Fe surface doping. Assuming Ising type easy axis order 
(moment || z) this leads to the opening of a mass gap 
at the Dirac point'' This may be derived from a 

RKKY mechanism including all spin componentd^. Fi- 
nally the last term is the Zeeman term with 75 = ^gsi^B 
where gs is the effective surface state g-factor which is 
in the rangd^ gs — 50. We note that under reflection 
fffi with respect to the horizontal disk center plane the 
Hamiltonian transforms as ah[H(k)] — iJ(— k). There- 
fore surface states on both sides of the disk contribute 
additively to M. We assume here that the disk is suffi- 
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ciently thick to prevent coupling between surface states 
on both sides. The evolution of surface states as func- 
tion of thickness was considered in detail in Ref. The 
warped Dirac cone dispersion is obtained from Eq. ^ as 



<=Dk 



-IsByf + {vFky+^sB^f (3) 



where we defined the warping function as 



(4) 



Contour plots of this dispersion in the {k^^ ky) plane that 
show the effect of warping , finite mass and Zeeman term 
are presented in Fig. [ij where intrinsic units m, B defined 
in the caption are used. The characteristics of the dis- 
persion will be discussed below. It is obvious that e^ik is 
the gapless isotropic Dirac cone dispersion shifted in the 
{kx^ ky)- plane to a new origin by the shift vector 



75 

VF 



{By, -Bx) 



IsB 

VF 



(sin 0, 



(5) 



In the following index 's' refers to shifted quantities. 
Defining the momentum q = k — kg with respect to the 
shifted orgin the dispersion in Eq. (|3| may be rewritten 
in the form 



= E, 
= E, 



Oq 



Oq 



T[4q2 + (TO + A$^)2 



(6) 



where the isotropy of the linear dispersive term (for 
small |q|) becomes obvious. Note, however that the 
parabolic and warping terms now depend manifestly on 
the magnetic field through the shift vector accord- 



ing 



to E^^ = (q +^ 

{{Qx + ksx)[iqx + ksxf 



,)V2m* and $^ 
-3{qy + k,yy]}. 




FIG. 3: Top: Torque oscillations as function of field angle 
(f} for warping A = 1. Full, dashed and dash-dotted lines 
correspond to rh = 0, 0.025, 0.05 respectively. The symmetry 
of fourfold oscillations ts{B, (f) + tt) — rg(_B, (j)) which appUes 
for m = is lost for nonzero rh. Bottom: Low field behaviour 
of torque at </> = 7r/8 . Same line conventions as on top. For 
nonzero rh asymmetry Ts{—B,(f>) ^ Ts{B,(j)) appears. Inset 
shows Dirac cones {rh = 0) with chemical potential position. 



where the surface states merge into the bulk. This rigid 
shift model is an approximation that does not take into 
account the details how the surface states merge into bulk 
at the cutoff and how this depends on the field induced 
shift of the origin. The parabolic term Uq leads to a con- 
stant contribution X^afAB) = — (7s/''^F"^*)'*S'(M)'^a0 ('^S 
— number of occupied surface states). It is proportional 
to the unit matrix and therefore does not contribute to 
the torque according to Eq. ([T]). Therefore, in this equa- 
tion we may replace Xap ~^ Xap since only the Dirac 
part Xa,3 may lead to a nonzero torque oscillation which 
is considered now. The leading part of the Dirac cone 
contribution per site due to Ud is given by 



For the torque effect we need to calculate the param- 
agnetic susceptibility associated with the above surface 
state spectrum which is defined by 
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(7) 



where Us is the contribution of surface states to the in- 
ternal energy per site. According to Eq. ^ it has two 
contributions 
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with = dH{p — EqT-) giving the occupation number 
of the surface state and fi denoting the chemical poten- 
tial which lies in the bulk gap. It should be stressed that 
the integration runs over the occupied and rigidly shifted 
Dirac states up to a cuttoff wave number Qc = fic/vp 



From now on we assume particle hole symmetry, i.e. ne- 
glect Eq in the Heaviside function 9h of the above equa- 
tion. Using J2t ■^^ff (m - TEq) = 9H{tJ'- Eq) -9H{tJ. + eq) 
one then obtains {Ac — A/N \s the surface unit cell area) 
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(10) 



Note that this integral extends over the occupied part 
(shaded area in the inset of Fig. [3]) for < 0. For fi > 
due to particle-hole symmetry the contributions from 
|eq.,.| < fi cancel and therefore Xq^(m) = Xq^(— holds. 
For 612X3 the cutoff qc is given by qc/qo — 0.1 with 
go denoting the zone boundary wave vector. The field 
derivatives are only nonzero due to the existence of the 
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warping term in Eq. ([6| and they are given by 



^dBj^dBp' 





"^^dBo^dBfs' 



(11) 



origin. For nonzero m this symmetry is lost (Fig. [TJl) 
for the dispersion and therefore also for ts{4') as shown 
later in Fig. [3j However, by definition the relation 

Ts{-B, + tt) = Ts{B, 0) stih holds. 



Evaluating the field derivatives of the warping function 



and defining $^ = 



and 



-) = ils/vp fqc^a/i (see Appendix 



.dB^dBfil - yia/'^-fJ 'ic-^qaji \°^^ x:^pi.^i^^i^|A| one fi- 
nally obtains the dimensionless susceptibility tensor of 
surface states as 



?/3 



q qafji 



de I dq _ 1 

[(j2 + (m + A$^)2]2 



(12) 



Here we introduced the dimensionless quantities 
q = q/qc, m = m//ic, A = {\/vF)ql = Ml/ tJ-c 
qp — (/i^ — m^)2 and B = ■ysB/fic- Obviously is 
only nonzero for finite warping parameter A. We define 
the amplitudes of the two torque contributions (Eq. ([T]) 
as Ai(0) = Ax^ = xfx - Xyy and A2(<^) = 2xSy. 



It is useful to consider first the low-field limit. For 



Xaf3 = with 



5 = the integral in Eq. (12) reduces to diagonal form 



q^dq 



[cp + rh?) 2 



(13) 



where 



Xo — /^c ^ls(lc/qo)^ ■ Therefore the zero-field 



susceptibility which vanishes identically for the isotropic 
Dirac spectrum becomes finite due to the effect of the 
warping. For rh = one has /xo = (9/5)7rA^(l — fi^) 
which is almost flat for fi < 0.4 and the drops to zero 
at /i = 1. However, as in the case of the parabolic 
Xaisi^) contribution x^^(B=0) is also proportional to 
the unit matrix. As a result the zero field torque ampli- 
tudes ^1,2(0) vanish meaning that the leading terms in 
^1,2(0) are of order ^ B^. Since in additon ^ B^ wc 
get rs((/)) ^ B^ for the low field behavior for the total 
torque. Explicitly, using Eqs. ( 1|12 ) the magnetic torque 
is then given by 



rs{<l>)^^^ = Ax^(0)sin2(/.-2xf,(</))cos2</) 



-VsLt^-'hsBfiK/ko)' 



(14) 



Here the overall amplitude Tg is determined by the vol- 
ume VsL of the surface layer, the Zeeman energy and the 
ratio of the Dirac cone cross section at /ic to the area of 
the surface BZ. 

Furthermore the torque amplitudes have to fulfill sym- 
metry requirements. In the massless case (m = 0) 
Eq. ([12]) and Eqs. ( |Al|A2[ ) lead to Ai,2(^ + 7r) = ^1,2(0) 
and consequently to ts{B, (f> + tt) = ts{B, (f>). This is in 
accordance with the twofold symmetry of the underly- 
ing surface state dispersion (Fig. [l|:) around the shifted 



For discussion of numerical results obtained from 
Eqs. ( |12|14D we give the typical magnitude of Dirac cone 
parameters for 612X3 in Table [ij Using values of fic and 
7s given there the reasonably accessible field range is 
< -B < 0.1 where the upper value corresponds to B = 
17.8 T. 



III. NUMERICAL RESULTS AND DISCUSSION 

The surface state dispersion given in Eq. ([3| is the 
fundamental property that determines the torque oscil- 
lations. It is shown in Fig.[I^-d (clockwise direction) by 
contour plots for different values of mass m and field B . 
There we used the intrinsic length scale b and energy scale 
E* related to the warping parametei^^. For {rh,B)— 
(0,0) we obtain the typical snowflake shape of the warped 
ungapped Dirac cones (a) that has been observed in 
photoemission"'. It exhibits both twofold and threefold 
rotational symmetry. The principal effect of the mass 
term is to remove these symmetries as seen in (b) where 
for clarity a large mass was chosen. The effect of the 
paramagnetic term due to an in-plane field and for zero 
mass is shown in (c): For finite B {(p ~ 37r/4) the center of 
the Dirac cone is now shifted from the origin of the BZ 
to the point {hk,j.,hky) = i3(sin — cos 0) — ^B{1,\) 
but the rotational symmetry around the shifted center 
is preserved at low energy. The direction of the shift 
depends on the direction of the field (chosen along a di- 
agonal here) . For finite mass and field both effects (sym- 
metry reduction and shift) are combined (d). 

The dispersion enters directly in the denominator of 
the susceptibility tensor in Eq. \12\ w hich determines the 
torque amplitudes yli^2 (</>)• In Fig. j2| they are first plotted 
as function of the warping parameter for a field angle (j) 
where both Ai_2 give non-zero contribution to ts{4>)- As 



TABLE I: Typical Dirac cone parameters for Bi?X.i(Ref. ITO]!. 
For the torque amplitudes the length and energy scales q^^ 
and fic are used. The intrinsic length and energy scales b and 
E* associated with the warping A are only used for Fig. [T] . 

Fermi velocity vp 2.55 eV A 

cutoff wave number qc 0.10 

cutoff energy /ic ~ vpqc 0.255 eV 

intrinisic energy scale E* — (vp/b) 0.25 eV 

intrinsic length scale b={X/vp)2 10.2 A 

warping parameter A = v^/E*^ 2.65 • 10^ eV A^ 

(dimensionless) X — \q"l/fj,c 1.04 

surface state g-factor gs 52 

field energy scale 7s = ^gsfJ-B 1-44 meV/T 
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FIG. 4: Top: Variation of = 7r/8 torque with chemical 
potential (fi = fi/fJ.c) for mass gap m — 0,0.05,0.1 (full, 
dashed, dash-dotted lines respectively). Due to particle-hole 
symmetry fs(/i) — fs(— /l). The minimum chemical poten- 
tial is = m. For \fl\ = 1 the Fermi level lies on the 
bottom of conduction or top of valence bands and the sur- 
face contribution fg vanishes. Bottom: Torque asymmetry 
Ats{B, <j>) = Ts{B, (j> + tt) — Ts{B, 0) as function of the mass 
gap for \ — 0.2,0.7,0.9 (full, dashed, dash-dotted lines re- 
spectively) 



can be seen directly from Eq. (12 1 the amplitudes are 
invariant under A — )• — A in the massless case (top) but 
not for finite m. This applies also for the total torque 
Ts (bottom). The Ai^2 amplitudes are nonmonotonic in 
A, however the total |ts(0)| still increases monotonically 
with |A|. For small warping and general <p '^si'P) ^ 
coming from the field derivatives of the warping function 
and for larger |A| fs{(f>) increases rapidly due to the effect 
of the warped dispersion in the denominator of Eq. ( 12 1. 



The value relevant for 612X3 is close to the intermediate 
value A = 1 which will be used in the following figures. 

The actual torque oscillations with field angle are 
shown in Fig. d (top) for i? = 0.1. They exhibit four- 
fold symmetry in the massless case (thick line) which is 
broken for finite mass. For the latter therefore a nonzero 
asymmetry Ats{B, </>) = ts{B, <P + tt) — ts{B, </>) appears. 
As can be seen in Fig. [s] (bottom) ts varies quadratically 
with B. This figure directly shows the asymmetry of the 
torque under 0—>(/)-|-7r or i3—>—_B for m ^ (broken 
lines) . 

In Fig. [4] we show the extremal torque amplitude as 
function of chemical potential. The torque oscillations 
become maximal when /t — !■ for half occupancy of the 
surface states. We remind that for /i > (chemical 
potential in the upper half cone in Figjsja) there is 
partial cancellation of contributions due to particle-hole 
symmetry and therefore Tsifl) = Ts{—fi)- The states 
close to the Dirac point show little warping, therefore 
when jl crosses that region the torque amplitude changes 
slowly. This is due to the warping function appearing 



in the numerator of Eq. (12). The strong variation of 
the torque for larger chemical potential in Fig. |4]jtop) 
would be the signature of its surface state origin in an 
experiment. This may possibly be achieved by an addi- 
tional gate voltage tuning on the disk as demonstrated 
in Ref. M 



IV. CONCLUSION AND OUTLOOK 

We have calculated the magnetic torque expected from 
helical surface states of a topological insulator when the 
magnetic field is rotated in the surface plane. For finite 
warping of the Dirac cones the paramagnetic moment in- 
duced in the surface states is non-coUinear with the field 
and leads to torque oscillations. The total amplitude in- 
creases with fourth power in the field. Furthermore the 
presence of a mass gap in the dispersion leads to a sym- 
metry reduction of the fourfold oscillations which should 
be detectable experimentally. The surface torque also de- 
pends considerably on the chemical potential position in 
the bulk band gap sufficiently away from the Dirac point, 
in contrast to the torque originating from the bulk states. 
Furthermore the ratio of the two contributions may be 
systematically changed by varying the layer numbers in 
thin film samples. These two effects should facilitate the 
separation of the surface torque from the constant bulk 
background. In this work we considered only constant 
and uniform fields. It seems possible to extend the the- 
ory to the case of dynamical and nonuniform torque using 
appropriate spin response function^^ll 



Appendix A 



Here we give the explicit form of dimensionless warp- 
ing function $| and its field derivatives , $|a/3 that 
appear in the numerator of the susceptibility tensor 



in Eq. (12). Using polar representation of momentum 
q = g(cos6',sin6') and field B = B {cos cf), sin (j)) coordi- 
nates we obtain for the first order derivatives 



{q cos 9 
\{q cos ( 



B sin (j) 
f B sin ( 



(Al) 

i{qsine ~ B cos (j)f] 



= i[q^ sin2e ~ 13^ sin2(j} - 2q]3 cos{e + (j))] 
= ^[f cos2e - B"^ cos2(t> - 2qB s\u{e + 4))] 



and likewise for the second order field derivatives 
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(A2) 



6 



In the case B=0 relevant for the zero-field susceptibility 



in Eq. ( 13 1 this reduces to 
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